Coulomb Charging Effects for Finite Channel Number 
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We consider quantum fluctuations of the charge on a small metallic grain caused by virtual 
electron tunneling to a nearby electrode. The average electron number and the effective charging 
energy are determined by means of perturbation theory in the tunneling Hamiltonian. In particular 
we discuss the dependence of charging effects on the number TV of tunneling channels. Earlier results 
for N ^> 1 are found to be approached rather rapidly with increasing N. 
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Single electron effects are well studied for metallic junctions in the region of weak tunneling In the case 
of a large tunneling conductance Gt comparable or even larger than the conductance quantum Gk — e 2 /h, it is 
necessary to go beyond lowest order perturbation theory in the tunneling Hamiltonian M - |hJ. In this case the 
theory involves contributions of order 1/iV where N is the number of transport channels. Since the junction area is 
typically much larger than X F , where is the Fermi wavelength, the number of tunneling channels N is often very 
large for metallic junctions, typically about 10 4 . Terms of order 1/N are thus dropped in most former approaches. In 
contrast to lithographically fabricated metallic junctions, in break junctions only a small number of channels may 
be available and terms proportional to 1/N cannot be neglected. Also for tunnel barriers jl2l in a two-dimensional 
electron gas there are typically only 2 spin degenerate transport channels contributing to the tunneling current. Thus, 
it is interesting to investigate how Coulomb blockade effects are modified by such 1/N corrections. Here we focus 
on systems where the dimensionless conductance per channel, go = Gt / Mtt JVGjr, remains small, that is to cases 
where the channels contributing to charge transfer are weakly transmitting. Further we restrict ourselves to the zero 
temperature case. This covers only partly the range of experimental interest but the results indicate how relevant 
1/N corrections are. 
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FIG. 1. a) Circuit diagram of a metallic grain (area within dashed line) coupled by a tunnel junction to the left electrode 
and capacitively to the right electrode, b) Electron-hole pair excitation in channel a of the tunnel junction. 



Specifically, we consider a small metallic grain in between two bulk electrodes biased by a voltage source U cx . The 
circuit diagram of the system is depicted in Fig. |l|a). Between the grain and the left bulk electrode electrons can 
tunnel while the right electrode couples purely capacitively. Thus no dc-current flows through the system. Due to the 
discreteness of the tunneling process, the excess charge of the grain can be shifted only by multiples of the elementary 
charge e, and we introduce an excess charge number n which characterizes the island charge q = —ne. As far as the 
tunneling conductance per channel is small compared to the conductance quantum Gk, electron tunneling can be 
described in terms of a tunneling Hamiltonian which will be treated as a perturbation. The unperturbed Hamiltonian 
for the Fermi liquids and the Coulomb energy reads 
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where n GX = GaU ex /e is a dimensionless voltage and 
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Ec ~2(C T + C G ) (2) 

is the charging energy needed to transfer one electron to the uncharged island at n cx = 0. eka and e qa are one- 
particle-energies for channel index a and longitudinal quantum number k and q, respectively. The operator a kcr (a,k a ) 
describes creation (annihilation) of a quasiparticle with quantum numbers ka, and a qr7 , a qa are defined likewise. With 
these operators the tunneling Hamiltonian may be written as 

Ht = Y, (>aL<vA + h.c.) (3) 

kqcr 

where t a is the tunneling matrix element and A shifts the island charge number n thereby changing the Coulomb 
energy. 

At zero temperature the system is described by the ground state energy £ and in view of eq. ([l]) the expectation 
value of the island charge may be written 

(n)=n ex -—- . (4) 
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To calculate £ we make use of Rayleigh-Schrodinger perturbation theory. The perturbation (|^) contains products of 
one creation and one annihilation operator, thus only even powers in t a contribute to £ , and it may be written 

£ = £^+£^+£^+0(4). (5) 

The zeroth order term £^ without tunneling is given by the minimum of the electrostatic energy and reads E c (n,Q — 
n C x) 2 where no is the integer closest to n cx . Hence, the averaged island charge in zeroth order perturbation theory 
is given by the well known Coulomb staircase. £ depends on n cx only via the electrostatic energy. It is thus an 
antisymmetric and quasiperiodic function of n cx which allows us to confine ourselves to < n cx < \. For the second 
order term we get formally 

£^ = {0\H T QH T \0) (6) 

with the auxiliary operator 

Q = ™ (7) 

In the same way the fourth order term reads 

£W = (0\H T QH T QH T QH T \0) ~ (0\H T Q 2 H T \0){0\H T QH T \0) (8) 

where terms with the energy denominator squared arise from the normalization of the ground state wave function. 
Inserting the tunneling Hamiltonian (|^) into the second order contribution (||) one gets 
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with the Coulomb energy differences 6E n = E c (n 2 — 2nn cx ). Both summands correspond to the virtual creation of 
an electron-hole pair with electron and hole sitting on different electrodes, cf. Fig. ^b). This can be represented in 
terms of Goldstone graphs depicted in Fig. |[ The contributions of second order correspond to the diagrams in Fig. 
||a) where the upper arc with an arrow to the right (left) electrode represents the creation of an electron-hole pair 
with the electron sitting on the right (left) and the hole sitting on the other electrode of the junction. The lower arc 
destroys this pair and we may omit the arrow since the process is uniquely determined by the upper one. 
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FIG. 2. Graphical representation of Rayleigh-Schrodinger perturbation theory, a) Graphs of second order, b), c), and d) 
Graphs of fourth order where graphs with inverted arrows are omitted. 

If we assume that the one particle energy is separable into a longitudinal and channel part, the summand in eq. (|]) 
is independent of the channel a apart from the factor t 2 . The sum over the transversal and spin quantum numbers 
leads to an overall factor N multiplied by the average t 2 = ^2 a t^/N of the tunneling matrix elements. We assume 
the bandwidth D to be large compared to E c and take the limit D/E c — > oo at the end of the calculation. To 
keep the finite bandwidth in intermediate formulas, we introduce an exponential cutoff and replace the sum over the 
longitudinal quantum numbers by 
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where p is the density of states at the Fermi level. Analogously, we introduce an integral with p' and D for the q sum. 
Formula <M) then takes the form 



£^ = -t 2 pp'N 
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We now define the dimensionless tunneling conductance per channel by go = t 2 pp' . Whereas the integral (|Tl]) is 
divergent for D — > oo, the average island charge to first order in go remains finite leading to 
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in accordance with previous results 




FIG. 3. Effective charging energy for channel number N = 5 as a function of the dimensionless tunneling conductance 
Gt/Gk in first and second order perturbation theory with and without 1/N corrections. 
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Along these lines the fourth order term m) is given in terms of double sums over kqa. These contributions can 
also be represented graphically as twofold virtual electron-hole pair creation and annihilation. We distinguish three 
groups of graphs depicted in Figs. ||b), c) and d). The first set of the graphs represents virtual electron hole pair 
creation and annihilation in two distinct channels a and a' which are only coupled by the energy denominator. Fig. 

represents graphs arising from the normalization of the wave function where the dash signs that the energy 
denominator is squared, cf. eq. (||). If there were no Coulomb interaction, i.e. E c = 0, the graphs b) and c) would 
cancel in accordance with the linked cluster theorem for uncorrelated fermions. Here we have to take these terms into 
account since the Coulomb energy correlates the electrons in the two electrodes. The channels of the graphs in b) 
and c) are not restricted and we get by summation over the transversal quantum numbers a factor N 2 . On the other 
hand, the graphs in Fig ^|d) describe processes within one channel because the electron created recombines not with 
the hole created at the same time but with a different one which has to have the same channel number. Thus we get 
by summation over the channels only a factor N. Therefore, we write the average island charge in the form 

f 1 4- In \ 

(n) = g N\n cx + (g N) 2 [c(n cx ) - c(-n ex )] + g 2 Q N[d(n cx ) - d(-n cx )} + O(g 3 ). (13) 

\ t -^cx / 

Previous theories for N 3> 1 introduce the dimensionless conductance g — goN of the junction and the terms 
proportional to g 2 N = g/N are dropped. 




FIG. 4. The average island charge number for dimensionless tunneling conductance Gt /Gk = 5 and channel number N = 5 
as a function of the dimensionless gate voltage n ex • 

The graphs in Fig. || b) and c) do not include 1/A-corrections and they lead to the known second order result for 
the island charge H 
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The graphs in Fig. || d) can also be integrated out analytically leading to a new contribution of order 1/N. We find 



d{u) = -In 3 



4 - 4m 



1 -2u 

[151n(2) + 161n(4-4u)]ln : 



26 q , 

ln 3 (l-2M) + 151n 3 (3-2i 

3 



l-2u 
3-2u 



[21n(l + 2u) - 71n(3 - 2u)\ ln 2 (l - 2u) 



10 In 
4 In 
6Li 3 



4 - Au 
1 - 2u 
l-2u 



10n 2 



+ 381n 2 (3-2?i) 



In 



3 - 2u 
l-2u 
3 - 2u 



3Li 2 
-8Li 3 



1 -2u 
3 - 2u 



3Li 2 



3 -2u 
A - Au 



1 - 2u 

3 - 2u 

-2Li 2 



8(1 - u) 

(3-2 U y 



3 - 2u 
A - Au 





2u\ 


8Li 3 ( 4 


- 2u\ 




Au) 




-Au) 



(15) 



where Li 2 (z) is the dilogarithm and 1^3(2:) the trilogarithm function fll3| . Within second order perturbation theory 
the result (|l3|) is valid for arbitrary channel numbers including N = 1 . 
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We now compare our result to earlier findings. One of the mostly frequently discussed quantities is the effective 
charging energy f| E* characterizing the effective strength of the Coulomb blockade effect. This quantity is defined 
by 



2dnL 



(16) 
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For small tunneling conductance, g — > 0, the effective charging energy approaches the bare charging energy E c 
whereas for strong electron tunneling E* vanishes. Our analytic expression leads to 

E*/E c = l-4g + 5.066 . . . g 2 - 7.167 . . . g 2 /N (17) 

where the constants are analytically known but too lengthy to present here. In Fig. |^ we show the normalized effective 
charging energy in first and second order perturbation theory without 1/N corrections compared to the complete 
second order result as a function of the dimensionless conductance Gt/Gk- We see that the 1/N corrections become 
more significant for larger tunneling conductance. 
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FIG. 5. Goldstone graph of order g§ which gives the leading asymptotic contribution for n cx — > A. 

In Fig. U the average island charge (n) in first and second order in go with and without the 1/N corrections is 
depicted for the case N = 5 and Gt/Gk — 5. We see that the 1/N corrections become significant especially for 



larger external voltages. From a comparison with Monte Carlo data |10| one estimates that results of second order 
perturbation theory are reliable for n cx up to 0.3. In the vicinity of the step, i.e., for n ex — > ±, finite order perturbation 
theory diverges and one has to sum diagrams of all order J^]. From our eq. ( |T5| ) one sees that the 1/iV-corrections 
become relevant near n ex = J even for large N since the qualitative behavior is changed. For n cx — > §, the leading 
N terms of second order show the logarithmic divergence || 2(goN) 2 In 2 6 where S = \ — n CK , whereas the 1/N 
corrections lead to — ^g^Nki 8 which eventually dominates the asymptotic behavior. For the special case N = 1, 
this divergence is in accordance with earlier findings by Matveev Q] . In general, for the finite N corrections of order 
<?g we get the leading asymptotic behavior 

(n) - g^iVln 2 ^ 1 S + 0{g k+1 ) (18) 

arising from the diagram depicted in Fig. |5| These terms dominate the asymptotics for n GX - ► hi even for large N . 

In summary, we have studied the influence of finite channel numbers on the Coulomb staircase. We found that 
for small voltages, n cx w 0, the 1/N corrections are significant up to N w 6 and become increasingly important as 
n cx = \ is approached. While near the step finite order perturbative results as derived here are not sufficient, the 
expression obtained for the effective charging energy is valid within an experimentally relevant range of parameters. 

This work was supported by the Deutsche Forschungsgemeinschaft (Bonn) and the Deutscher Akademischer Aus- 
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